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Abstract 

The purpose of this paper is to investigate the spectral nature of the Neumann-Poincare 
operator on the intersecting disks, which is a domain with the Lipschitz boundary. The 
complete spectral resolution of the operator is derived, which shows in particular that it 
admits only the absolutely continuous spectrum, no singularly continuous spectrum and 
no pure point spectrum. We then quantitatively analyze using the spectral resolution the 
plasmon resonance at the absolutely continuous spectrum. 
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1 Introduction 

The Neumann-Poincare (NP) operator is a boundary integral operator naturally arising when 
solving classical Dirichlet or Neumann boundary value problems using layer potentials. This 
operator is not self-adjoint with respect to the usual inner product unless the domain is a 
disk or a ball. However, it is recently found that by introducing a new inner product the NP 
operator can be symmetrized and admits a spectral resolution. If the boundary of the domain is 
smooth, then the corresponding NP operator is compact, and hence has only the point spectrum 
accumulating to 0. But, if the boundary of the domain is Lipschitz, then the NP operator is 
not compact (it is a singular integral operator), and it is quite interesting to investigate the 
nature of the spectrum of the NP operator, in particular, whether it admits both point and 
continuous spectrum. In this paper we derive a complete spectral resolution of the NP operator 
on the intersecting disks which is a Lipschitz domain, from which we are able to show that the 
NP operator on the intersecting disks has only absolutely continuous spectrum. As far as we 
are aware of, this is the hrst example of a domain with the Lipschitz boundary for which the 
complete spectrum of the NP operator is derived. 

The study of the NP operator goes back to Poincare as the name of the operator alludes. 
It has been a central object in the theory of singular integral operators developed during the 
last century (see m)- Recently there is growing interest in the spectral property of the NP 
operator in relation to plasmonics. The plasmonic structure consists of an inclusion of the 
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negative dielectric constant (with dissipation) embedded in the matrix of the positive dielectric 
constant. The eigenvalue of the NP operator is related to the plasmon resonance (see for example 
0)- The negativity of the dielectric constant of the inclusion makes it necessary to look into 
the spectrum of the corresponding NP operator. If the dielectric constant is positive, which 
corresponds to the ellipticity of the problem, only the resolvent of the NP operator matters. In 
this paper we derive quantitatively precise estimates of the plasmon resonance at the continuous 
spectrum of the NP operator on intersecting disks in the quasi-static limit. Particularly we 
obtain precise rates of resonance as the dissipation parameter tends to zero. As far as we are 
aware of, this is the hrst such a result. We mention that the plasmon resonance (with a non-zero 
frequency) at the continuous spectrum on the intersecting disks was also studied in [15]. But 
there quantitative estimate of resonance is missing. 

As mentioned before the NP operator is not self-adjoint with respect to the usual inner 
product, unless the domain is a disc or a ball (see m)- However, it is recently found in m (see 
also [9|) that NP operator can be symmetrized using Plemelj’s symmetrization principle. There 
are a few domains with smooth boundaries on which the spectrum of the NP operator is known; 
a complete spectrum is known on disks, ellipses, and balls, and a few eigenvalues on ellipsoids 
(see [9] and references therein). The spectrum of the NP operator on two disjoint disks has 
been computed recently, and applied to analysis of high concentration of the stress in between 
closely located disks laEiiiT]. The spectral property of the NP operator is barely known. But, 
the bounds on the essential spectrum of the NP operator on the curvilinear polygonal domains 
have been derived in m- The intersecting disks is among domains considered in the paper, 
and interestingly the result of this paper reveals that the bounds there are optimal for the 
intersecting disks. Recently, the spectral theory of the NP operator finds new applications. In 
[T] it was efficiently used for analysis of the cloaking by anomalous localized resonance on coated 
plasmonic structure. More recently, uniformity (with respect to the elliptic constants) of elliptic 
estimates and boundary perturbation formula has been proved using the spectral property of 
the NP operator Unj. It is also applied to analysis of plasmon resonance on smooth domains [1]. 

The finding of this paper shows that the NP operator on intersecting disks does not admits 
a point spectrum on ' (dQ). We suspect that this is rather exceptional and domains with 
corners, especially 4 or more corners like squares, admit point spectrums. This is an interesting 
question to investigate. In this regards, it is also interesting to see whether all the non-smooth 
Lipschitz domains admit non-empty continuous spectrum. It is worth mentioning that the NP 
operator on space on any domain Q (with possibly multiple connected components) 

with the Lipschitz boundary has 1/2 as an eigenvalue and its multiplicity is the same as the 
number of connected boundary components [2]. 

This paper is organized as follows. In the next section we introduce the NP operator and 
the bipolar coordinate system which is an essential ingredient of this paper. In section [3| we 
explicitly compute the single layer potential and the NP operator on the intersecting disks in 
terms of the bipolar coordinates. In section H] we derive a resolution of the identity and the 
spectral resolution of the NP operator. The last two sections are to derive estimates for the 
plasmon resonance. 

2 Preliminaries 

Let us fix some notation first. Let 11 be a bounded domain in with the Lipschitz boundary. 
We denote by the dual space of where the latter denotes the usual Sobolev 

space. We use the notation (•, •) for either the L?' inner product or the duality pairing of 
and on 912. We denote by || • ||s the Sobolev norm on for s = —1/2 or 1/2. Let 

be the space of V’ £ satisfying {ip, 1) = 0. 
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2.1 The NP operator and symmetrization 

For X = {xi,X 2 ) / 0, let r(a;) be the fundamental solution to the Laplacian, i.e., 

r(x) = — In |x|. (2.1) 

27r 

A classical way of solving the Neumann boundary value problem on P is to use the single layer 
potential Sqq [^p] which is defined by 

SavLVp\{x) ■■= ( T{x-y)p{y)da[y) , x G (2.2) 

Jan 

where p G Lp'^dVl) or It is well known (see, for example, laiz]) that satisfies 

the jump relation 

■^Saniv’ll^ix) = (^±^1 + ICgQ^[(p]{x), a.e. X G , (2.3) 

where the operator ICqq is defined by 

^ani^Kx) = ^ • (2-'^) 

Here ± indicates the limits (to 9H) from outside and inside of H, respectively, and i^x denotes the 
outward unit normal vector to dQ at x G dP. The operator ICqq is called the Neumann-Poincare 
(NP) operator associated with the domain P. Sometimes it is called the adjoint NP operator to 
distinguish it from the NP operator ICan {^an is the adjoint of K,*q^)- 

Even though 1C*qq^ is not self-adjoint on the usnal L^-space, it can be symmetrized by intro¬ 
ducing a new inner product. Recently it is proved in |14] that can be symmetrized nsing 
Plemelj’s symmetrization principle (also known as Calderon’s identity) 

San^an = ^anSan- (2.5) 

In fact, if we define 

■=-{p,San[i>]) ( 2 . 6 ) 

for ^ Hq ' (dP), then it is an inner product on Hq ' (dP) and K,*q^ is self-adjoint with 
respect to this inner prodnct. 

Let Ti* = 'R*(dP) be the space Hq ' (dP) equipped with the inner prodnct ( , )'^* and 
denote the norm associated with ( , )-^* by || • Wu*- It is proved in [TO] that || • ||-^* is equivalent 
to norm, i.e., there are constants Ci and C 2 such that 

Cl 119^11-1/2 < WfWu* < C'2||(/?||_i/2- (2.7) 

for all p G Hq 

Since is self-adjoint on %*, its spectrnm ct^ICqq^) is real and consists of point and con- 
tinnons spectra. Moreover, by the spectral resolution theorem there is a family of projection 
operators St on T-L* (called a resolution of the identity) snch that 

/■ 1/2 

IC*g^= / tdSf (2.8) 

J- 1/2 

See [21] • We emphasize that the spectrum of ICqq onTJ* lies in (—^,1) (see [12] )■ 
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2.2 Intersecting disks and bipolar coordinates 

Let Bi = Ba{ci) and B 2 = Ba{c 2 ) be two disks of the same radii a centered at ci and C 2 , 
respectively, which are intersecting, namely, |ci — C 2 I < 2a. The domain Q we consider in this 
paper is defined to be 

n:=BiUB2. (2.9) 

(We may consider Bi n B 2 instead.) Note that dQ is Lipschitz continuous with corners at the 
intersecting points. We assume that ci is on the a; 2 -axis and ci = —C 2 so that the intersecting 
points lie on the xi-axis and they are symmetric with respect to the origin. Let 26q {9q > 0) 
be the angle between two circles at the intersecting point (see Fig. 1 ^. One can easily see that 
the intersecting points are given by iasin^o- 



-5 0 5 



Figure 2.1: The left figure is the intersecting disks 11. Level coordinate curves of C (when 
a = \/7) are illustrated in ^;-plane in the right figure. 


For a given positive number a the bilinear transformation 

A(.2) = (2.10) 

z — a 

is a conformal map from the Riemann sphere onto itself, and maps the interval [—a, a\ onto the 
negative real axis. So, 

C =Log A(z) (2-11) 

is well defined for z £ C \ [—a, a], where Log is the logarithm with the principal branch. The 
coordinate C, or (^, 6) is the bipolar coordinates we will use. The level coordinate curves for C, 
are illustrated in Fig. l2.1l Here and afterwards, we take 


a = a sin 6 * 0 , 


( 2 . 12 ) 


so that (a,0) and (—a, 0) are two intersecting points of the circles. Note that LogA maps 
C \ [—a, a] onto {A = ^ + i6, —00 < .^ < 00, —vr < 6 < n}. If we let F_|_ and r_ be the lower and 
upper halves of dQ, respectively, one can see that LogA maps F_|_ onto {9 = and F_ onto 
{9 = —0o}) and outside of 11 onto {A = ^ + i9, —00 < ^ < 00, —9o < 9 < 0o} (see Fig. 12.21) . Note 
that we denote the lower half by F+ because it is mapped to {9 = +0o}- 
Let <I> be the inverse of LogA, he.. 


Z = Xl+ iX2 


<I>(C) :=a 


+ 1 
el — 1 


(2.13) 
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Figure 2.2: The exterior of the intersecting disks (the gray region in the left figure) is the image 
of the infinite strip (the gray region in the right figure) via the conformal map z = xi + ix 2 = 


We also regard ‘k as a mapping from into by identifying x = (xi,X 2 ) S with 2 = 
xi + ix 2 G C. Here we recall a few facts to be used later. We define the scale factor h by 


1 cosh,^ — cos 0 

(2.14) 

Since <1> is a conformal transformation, h is the Jacobian of the coordinate change. In particular, 
one can see that the outward normal derivative on 90 is transformed in the following way: 

8„ = =f m 

(2.15) 

and the line element da on 90 is transformed by 


HLGq) 

(2.16) 

Define 

LogA(^;) =: Ti(z) + i^ 2 {z)- 

(2.17) 

Then we have 

f$(Ti(x),6»o) = X ifxGr+, 

|$(Ti(x), — 0o) = a: if x G r_. 

(2.18) 


3 Some explicit computations 

3.1 Single layer potential on 

The purpose of this section is to derive an explicit formula for the single layer potential on 90 
using the bipolar coordinates. For (p S T-L* let 

u{x) := Sdn[(f]{x), x G M^. 
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Then u is the solution to the following transmission problem: 


Au = 0 


u\ + 
du 
du 


= u\- 
du 
+ dv 


= ^ 


u{x) = 0{\x 




in n U (M^ \ 0), 
on dn, 

on dn, 
as \x\ oo. 


(3.1) 


It is worth mentioning that the last condition follows since {ip, 1) = 0. 

Let us define 

u{^,0) := {uo^){^,6), 

where <I> is the conformal mapping dehned by (j2.13l) . Since is conformal, u satisfies 

Au{^,e) = 0 for (^,0) e Mx((-7r,-6'o) U (-6»o,6»o) U (6»o,7r)), (3.2) 


where A = 9^/(9^^ + jdO^. The continuity of the potential in (13.ip is equivalent to 


u\+ = u\- on {6 = ±0o}- 

In view of (I2.15p . the continuity of the flux is equivalent to 


du 

dd 


and 


du 

dd 


du 

” de 

du 

de 


-1 


- h{;9o) 

1 


(p+ on {e = 00 } 


ip_ on {9 = -0o}, 


(3.3) 

(3.4) 

(3.5) 


- h{;9o) 

where 

<^+(0 = ^)(C>^o) and (p-{0 = {ip o -do). 

Since |(xi,X 2 )| —^ oo if and only if |(C,0)| ^ 0, it follows from the last condition in (13.ip that 
u{^,9) —)• 0 as {i,9) —)• (0,0), or simply 


h(0,0)=0. (3.6) 

We now impose boundary conditions on {9 = Tvr}. Since lines {9 = ±7r} correspond to the line 
segment {—a,a) in (xi, X 2 )-plane, we have 

/ 77 / fjqi 

u{^,Tr) = u{C,-7r), for all ^GM, (3.7) 

which implies that u{^, 9) is periodic in 9 with the period 27r. 

It is straightforward to find the solution u to (j3.2p satisfying (I3.3p - (j3.7l) . Denote by T the 
Fourier transformation in .^-variable, namely, 

V J —QQ 

If we apply T to the both sides of the equation ()3.2I) and denote 

u(s, 0) = T'['u(-, 0)](s), s G M, 


then we have 


d^v 2 n 
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We solve this ordinary differential equation with conditions (I3.3I1 - (I3.7I) . In doing so, it is conve¬ 
nient to consider separately the cases when y? is odd and even with respect to xi-axis. 

Suppose that ip is odd with respect to xi-axis, namely 

p{xi,X2) = -p{xi, -X 2 ). 

Then, one can see easily that 

ip+ = -ip_, (3.8) 

and the solution is given by 

60 < 9 < TT, 

-9o<e< do, (3-9) 

—vr < 6 < —60, 

{s). 

If p is even with respect to xi-axis, then we have 

ip+ = (p_, (3.10) 


v{s, 9) = < 


a(s) sinh s{9 — vr), 
smh s9q 

a{s) sinh s{9 + vr). 


where 


, , sinhs^o 
os = ——- J- 


ssinhsvr L/i(-,0o)- 




and the solution is given by 


9) = (const.) -|- < 


h{s) cosh s{9 — vr), 
cosh s9o 

h{s) cosh s{9 + vr). 


9q < 9 < -K, 

-9o<9< 9o, 
—vr < 9 < —9o, 


(3.11) 


where 


coshs6»o ^ 
b[s) = -r-:- J- 




is). 


ssinhsTT \-h{-,9o)- 
In general case, we decompose p into the odd and even parts, namely, 

p = p° + p^^ 


(3.12) 


where 

V?°(xi,X 2 ) := ^{p{xi,X2) - p{xi,-X2)) and p‘'{xi,X2) := ^{p{xi,X2) + p{xi, -X2)), 

and then obtain the solution by superposing solutions corresponding to the odd and even parts, 
respectively. By taking the inverse Fourier transform we obtain u{^, 9), namely, 

u(C,0) = n($(C,0)) = J^-'K,0)](O. (3.13) 


We emphasize that 


In fact, since p G Ti*, it holds that 



p^da = 0 . 


Since p° is odd with respect to the 
p^ is even, we have (|3.14l) . 


/r/ 


p°da + / p^^da = 0 . 


xi-axis, we have Jp p°da 


(3.14) 


0, and hence fj. p^da = 0. Since 
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3.2 Single layer potential on dfl 

In particular, we obtain from (13.131) the following formula on 911: 


where 






— sinh sOq sinh s(7r — 6 * 0 ) 


s sinh STT 

K-M- 

— cosh s9q cosh s(7r — 6 * 0 ) 

\ 

s sinh STT 

Lm-,0o) 


(3.15) 


It is worth reminding here that 
and likewise for 

Above formula suggest to introduce an operator C as follows: for ip ^T-L* 


cwm ■■= 


c+[^]{C) 

c-m). 


where 

■= HKM ■ 

Note that C~^ is given by 


c-^ 

Qi 

|/i(Ti(x), 6(0)52 (^i(®)), xGT-. 

(3.16) 


[52 J 

We also define a multiplication operator by 




Mh[f] :=MC,0o)/(O- 


Let 


, , sinh S 0 O sinh s( 7 r — 00 ) 

(3.17) 



1^1 « — ■ V. 

s sinh 57 r 



, , cosh S 0 O cosh s( 7 r — 0 o) 

P2 5 ) := 

s sinh 57 r 

(3.18) 


Observe that pi is a smooth function on and p 2 is smooth except at 0 where the following 
holds: 

lim |spp 2 ('S) = 1. (3.19) 

s ->-±0 

Both of them are positive and satisfy 


lim \s\pj{s) =-, j = l,2. 

s—>-±oo 2 


(3.20) 


Then, since ¥?+(0 = —^-{0 (f'piO = P-{0 j see from (13.1511 that 

TMhC+Son[^]is) = -pi{s)J'C+[p°]{s) - p 2 {s)J'C+[p^]{s), 
TMhC-SQn[p]{s) = -pi(s)T’C_[v9°](s) - p2{s)TC-[p^]{s). 





















Let 

Then we have 
in other words, 

where 

So, we obtain 


/i(s)’ 


'fC+[p%s) 

Ms) 


FC+[p%s)_ 


f(s) = 

fi + f 2 = [ip] and - /i + /2 = FC- [ip ], 

p = 2C-^F-^Ar^ 


(3.21) 


Note that 


where 


A = 


2FM.hCSoQ{AFCy 


-Pifi -P2/2 
Pifi -P2/2 


1 -1 

1 1 


(3.22) 

(3.23) 


1 

7 i' 


1 

1 

to 


A 


> 

1 

to 


= -2A-^P(s)f(s) 


P = 


Pi 0 
0 P2_ 

If we define the transformed single layer potential S by 


it follows that 


(3.24) 

(3.25) 


§ := kFMhCSan{AFC)-\ 

S[f](s) = -P(s)f(s). 

3.3 Inner product 

Observe that if we make changes of variables x = <1>(^, ±0o), then we see from (12.161) that 

[ p{x)'tp{x)da{x) = [ C+[p]{C)MhC+[tp]{C)dC+ [ C-[p]{C)MhC-[tp]{C)dC- 

J do. Jm. Jm. 

It then follows from Parseval’s identity that 

f p{x)y{x)da{x) = ( FC+[p\{s)FMhC+[ip]{s)ds + [ FC^[p]{s)FMhC-[y]{s)ds 

Jao Jn. JK 

= [ FC[p]{s) ■ FMhCms)ds. 

Jr 

Since A“^ = ^A^, we have 

f p{x)'y{x)da{x) = - f AFC[p]{s) ■ AFJAhC['tp]{s)ds. 

Jao 2 

f = (/i,/ 2 )^ := AFC[p\ and g = ( 51 , 52 )^ := AFC[y]. (3.26) 


Let 
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Then, we have from (13.251) 




-{(p,Sdn['il;]) 


-- / ATCMs) ■ ATMhCS9nms)ds 
^ Jr 

-If i{s) ■S[g]{s)ds 
^ Jr 

I [ i{s) ■ P{s)g{s)ds. 

^ Jr 


Let 


1 

2 


f(s) • P{s)g{s) ds 


1 

2 


Pi{s)h{s)gi{s) +P2{s)f2{s)g2{s) 


ds, 


(3.27) 


and let || • \\w* be the norm induced by this inner product. Let W* be the space of f such that 
||f||w* < oo. In fact, W* = L^(M,pi) x L^(M,p 2 )) where L‘^(R,pj) is the space weighted by 
Pj. It is worth mentioning that because of (I3.20p L^(M,pi) is the Fourier transform of the space 
77“^/^(M). On the other hand, because of (|3.19l) . if 5 G L^(M,p 2 )) 5(0) = 0 if g is continuous at 
0. So, L^(M,P 2 ) is the Fourier transform of 77g ^^^(M). This is natural in view of (I3.14p . 

Define 

U := AFC. (3.28) 

We have shown that 1/ is a unitary transformation from P* onto W*, i.e.. 


= {U[(p],U[lp])w*- 


(3.29) 


3.4 Neumann-Poincare operator 

We now compute the NP operator and show that 

IC*Q^[p]m,±eo)) = h{C,eo)F~^[±B° - (3.30) 


where 


B%s) := r]{s)F 
B^s) := p{s)F 




h{;eo)\ 


(•s), 

(s). 


Here, g is defined by 


1 sinh s( 7 r — 29q) 


(3.31) 


We only show (13.301) for 6 = 6q. Let u = 590 ( 99 ] and u{^,9) = 59 o[( 93 ]($(^,0)) as before. We 
see from (12.3p and p2.15l) that 


/CS^[99](<i>(e,0o)) = 


1 / du 


2\dv 


+ 


du 


du 

h{^,0o) 

de 


imoo)) 


(C,0o) + |^|_(C,0o) 


(3.32) 


Recall that the region \ D outside the intersecting disks corresponds to the region {|5| < 0o}- 
So the subscript + in indicates the limit from the region {|5| < 0o}- Similarly, the subscript 
— indicates the limit from the region {Oq < \9\ < vr}. 
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Let V = v° + := F[u\. One can easily see from (|3.9I) and (|3.1ip that 


dv° sinhs(0o — tt) coshs^o -r- 

{s,9o) = -—- F 

+ sinh svr 


do 

dv° 

'W 

dv^ 

dv^ 




Lh(.,0o)J 


, cosh s (00 — tt) sinh S0O ^ 

- ’ ^ sinhsvr L/i(-,0o)- 

, cosh s (00 - vr) sinh S0O -r-l 'f+ 

(s,0o) =-—- F 

sinhsTT 

. sinh s (00 - tt) cosh S0O ' 

(s,0o) =-—- F 

sinnsTT 


(s), 

(s), 


(s). 


Since 


we obtain 


sinh s(0o — vr) cosh s0o + cosh s(0o — tt) sinh s0o = — sinh s(7r — 20o), 


dv 

de 


/ ^ N dv , ^ , sinh s(7r — 20o) , ^ 

+ 7 ^ {s,eo) = - 

o9 - sinhsTT 




Vh{-,9, 


oj-i 


(s)-^ 




Lh(-,0o)J 


(^) 


So, we obtain (13.301) for 0 = 0o. 

Dehne f by (13.211) . Then, (13.301) can be written as 


2FClCl^{kFC)-^[i] = r^ 


/i - /2 ' 

-h - f2 


= 277 A 


-1 


fi ■ 
-/2 


Let 


K* := U}C*9^U-^ = AFC}C*a^{AFC)-\ 

Then using (13.22^ . (I3.30p can be rewritten as 


K*[f](s)=r?(s) 


/i(s) 

-/2(s) 


on VT* 


(3.33) 

(3.34) 

(3.35) 


Note that rj is monotonically decreasing for s > 0, tends to 0 as s —)• 00 , and satishes 

20o(7r-0o) 2 


r/(s) = b 


1 - 


-s^ + 0{s^) 


(3.36) 


where 


h ■= ?7(0) = ^ 


(3.37) 


We can infer from (j3.35p that [—b,h\ is the spectrum of K* (and hence of Fqq), and it is a 
continuous spectrum. We will investigate the nature of the spectrum in the following section. 


4 Spectral resolution of NP-operator 

4.1 Resolution of the identity 

Let Xs be the characteristic function on (—00, s] for each s G M. Let 


Rf = 


and 112 f = 
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We define a pair of operators V] and on W* by 

Pi[f]=y,n/, j = l,2. 

One can see that Vi is an orthogonal projection, namely, Vi is self-adjoint and {Vi)"^ = Vi. We 
also have 

VlVt = VtVl = 0 for all s,t G M. (4.1) 

One can also see that the following holds: 


lun ^^If] = n,t, 
lim P||f]= 0 , 

oo 

limP/[f]=P^[f], 

t^s 

[vivi[i] = viM 


(4.2) 


where t A s denotes the minimum of t and s, and the convergence is in W*. As a consequence, 
we have 

lim (Pi[f]+P2[f]) = f. 


We see from (|3.27l) that 


{f,Vi[f])w* = ^ J Pj\fj\‘^ds, j 


= 1 , 2 . 


(4.3) 


(4.4) 


It then follows that 


(f,f)w* = ;^ [ [pil/iP+P2I/2P] ds. 
^ Jr 


roo j r°° d 

By the change of variables t = r]{s) for s > 0 and t = r](—s) for s < 0, we obtain 

f'b 

i(f, (pL-i.o -Pi-wn) m\ 

which can be written as 


(4.5) 


/Z (pVh.) -Pl-w) m 


r * = /Zii’ (A,-.,,) - A-«) +^) m 


w* 


where X is the identity operator on W*. If we make change of variables —t = r]{s) for s > 0 and 
—t = r]{—s) for s < 0, then we have 


/I *=//('’ (a'-w - a,-h<)) m 

So we obtain from (14.5p that 

(pVh.) - A-(.)+^) m)„,. 


w* 


(4.6) 
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We now define a 

family of projection operators {Et}jg[_b,6] on W* as follows: 



- K-pt) + tei0,b], 

(4.7) 



and 

En := lim E^, 

(4.8) 


4-:-0+ 

in other words, 

Eo[f] = lim Et[f] 



t-:-0+ 


for f G W*. In view 

of the first two identities in (14.2p. we have 



Eo[f] = -nif + f = n 2 f. 

(4.9) 


Now one can easily see from (j4.1h . (I4.2h . and (14.8h that the following holds: 

EtEs =EtAs, 

< Eh[i] = f, E_b[f] = 0, (4.10) 

lim Et[i] = E,[f] 

^ t—>-s + 

for all s, t G [—b, b]. 

We obtain the following proposition. 

Proposition 4.1 The family is a resolution of the identity on W*, namely, it holds 

that ^ 

= j ^d{g,Et[i])w* (4.11) 

for all g, f G W*. Moreover, we have 


limEi[f] =E,[f] 


(4.12) 


for all s G [—6, 6]. 

Proof. That {Et}jg[_j b] is a resolution of the identity on W* is an immediate consequence of 
(I4.10p . To prove (I4.12p . it suffices to consider the case when s = 0 because of the third identity 
in (j4.2p . But one can see from the definition ()4.7p that 

lim EJf] = lim Vt[{] - lim Vf[{] = Bsf. 

i—>- 0 — t—>00 t—^—co 

So by (14.op we have ()4.12l) for s = 0, and the proof is complete. □ 

4.2 Spectral resolution of the NP-operator 

We now derive the spectral resolution of IfC*. 

Theorem 4.2 Let {Ej} be the resolution of the identity defined in (14.7h 
following spectral resolution o/K.* on W*: 

K* = y tdEt. 


Then we have the 

(4.13) 
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In other words, it holds that 


(g,K*[f])w^* = r td{g,Et[I])w* (4.14) 

J-b 

for all g, f G W*. 

Proof. It suffices to show (j4.14l) when g = f by the polarization identity. We obtain from (I3.27P 
and (j3.35p 

(f,K*[f])vi/* = l [pil/ip-P 2 |/ 2 n ds. 

^ JR 

It then follows from (|4.4I) and (|3.31l) that 

/ CXD r I-OO J 

By applying the same changes of variables as before, t = r/(|s|) or t = —r/(|s|), we obtain 

pb rO 

= / td{f,Eti)w* + / td{i,Etf)w*- 

Jo J-b 

Since is continuous at 0, we obtain (|4.14l) when g = f. This completes the proof. □ 

4.3 Spectrum of the NP operator on the intersecting disks 

Recall that the point spectrum s is characterized by the condition (see |21j i 

lim Ef 7 ^ Eg. 
t — ys — 


Since E^ is continuous for all s, namely, (I4.12p holds, we infer that there is only continuous 
spectrum and no point spectrum of K*. In this section we prove that there is only absolutely 
continuous spectrum. For that we derive the following lemma. 


Lemma 4.3 For all f, g G W*, it holds that 


^(f,Et[g])w* 


Pi(g) 

2|r/'(s)| 

P2{s) 

.2|r/'(s)| 


fi{s)gi{s) + fi{-s)gi{-s)^, 

f2{s)g2{s) + f2{-s)g2i-s)\ , 


t G (0, b], 
t G [—b, 0), 


where s = g ^{t) if t > 0, s 


g ^{—t) if t < 0 and the equality holds in L^{[—b,b]). 


(4.15) 


Proof. By the polarization identity it suffices to consider the case when f = g. By the definition 
(|4.7p of the resolution of identity {E^}, we have 




- (f,^Vi(_i)[f])vu*, t G [-6,0). 
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So, if t > 0, we have from (14.411 that 


1 rv hi) 

(f,Ei[f])H^*=(f,f)^,-- / Pil/il 


' ds. 


Since pi is even, we obtain 
d 


dt 




Since rj' is negative, we obtain (I4.15P for t > 0. The case when t < 0 can be proved similarly. 

□ 

One can easily see that is integrable on [—6,6]. In fact, since f G W*, we have 


/ 


b 

-b 




1 

^*= 2 ]^ PiCs) (l/i('S)P + l/i(-s)n ds 

1 r°° 

^2 Jo (l/2(s)P + l/2(-s)n ds 


< oo. 


Since K* is unitarily eqnivalent to fCg^, we obtain the following theorem. 


Theorem 4.4 Let daciK^gfi), (Tsc{KLqq) o.nd asc{^*QQ) be the absolutely continuous, singularly 
continuous, and pure point spectrum of ICqq on Li*, respectively. Then we have 

O'aciK^an) — [~b,b], CTsciJCgQ) = 0, 0 'pp{ICgQ) = 0. 


A few remarks on Theorem 14.41 are in order. If 0o = i>'/2 or fl is a disk, then 6 = 0 which 
is in accordance with the fact that JCg^ = 0 on LL* when 0 is a disk. It is also interesting to 
observe what happens when 0o ^ 0; or when the intersecting disks approaches to the tonching 
disks. See Fig. 14.11 In this case the spectral bound tends to 1/2. This fact was also observed in 
the other direction in It is proved that if two separating disks approach to each other, 

then there are eigenvalues approaching to 1/2. It should be mentioned that the NP operator for 
the touching disk is a hyper-singular operator due to the cusp, and is not well-defined on 

The analysis of this paper goes through with slight modification for the case when Oq > tt/2 
(see Fig. 14.ip . In this case the spectral bounds is given by 


1 _ TT - 6>o _ ^ _ 1 

2 TT TT 2 


(4.16) 


Recently a bound on the essential spectrum of the NP operator on curvilinear polygonal domains 
has been obtained in [I9]. It is quite interesting to observe that the bound in (I4.16P is exactly 
the one obtained there for the intersecting disks. 


5 Plasmonic resonance on the intersecting disks 

In this section we consider plasmonic resonance on the intersecting disks. Suppose that the 
intersecting disks has dielectric constant Cc + id (6 is the dissipation) and the matrix \ has 
the dielectric constant Cm- So, the distribution of the dielectric constant is given by 

e = {ec + id)x{^) + emX(E^ \ 
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Figure 4.1: Disks of various intersecting angles 


where x denotes the indicator function. After normalization we assume that = 1- To 
investigate the resonance on the plasmonic structure, we consider the following problem: 


V • eVu = / in 

u{x) = 0(|x|“^) as |x| ^ oo, 

where / is a source function compactly supported in \ D and satisfying 

/ = o. 


/ 


(5.1) 


(5.2) 


A typical such source functions are polarized dipoles, namely, f{x) = a ■ V6z{x) for some 
2 ; G \ D, where a is a constant vector and 6 z is the Dirac mass. It is worth mentioning that 
(j5.2p is necessary for a solution to (j5.ip to exist. For a given Cc resonance is characterized by 
the fact 

l|Vu5||i2(t ^)00 as 50. 

Let 

Cc “1“ 1 “k iS 


2 {ec-l) + 2i6' 

Then, the solution us to dSI]) can be represented as 

us{x) = q{x)+SdQ[ips]{x), X G R^. 
where q is the Newtonian potential of /, i.e., 

q{x)= / T{x - y)f{y)dy, 

jRd. 

and the potential ips G Ti*{did) is the solution to the integral equation 

{XI - )C9Q)[ips] = dyq on 9D. 

See m for the proof of this representation. 

Let 

£t := U-^EtU 


(5.3) 

(5.4) 


(5.5) 


(5.6) 


(5.7) 


where Ef is the resolution identity on W* (Proposition 14.11) and U is the unitary operator in 
(j3.28p . Then we see from Theorem 14.21 that 


^dn — / tdSf 
l-b 


(5.8) 
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We then obtain from (15.61) 


= 


f 


1 


■d£t[dyq]. 


'-6 A - t 

Since ||V(;||^ 2 (q) < oo, we infer from (j5.5l) that ()5.3p is equivalent to 

||V5ao[</?5]||i,2(0) = W'^ius - q)\\L2(n) oo as <5 0. 

We see from (j2.3[) that 


r g _ 

||V5aQ[(^5]||^2(Q) = ■^Ssn[‘Ps] _5ao[‘/?<5] da 


Using ()5.9p we see that 

Ws, d{^s,£t[^5])'H* 

So, we infer that 

CiWifsWn* < ||V5ao[v95]||2,2(f7) <C2\\ips\\H* 
for some positive constants Ui and C 2 , and (I5.10p is equivalent to 

W^sWw -^00 as 5 0. 

It follows from (15.91) that 


/ 




d{dyq,£t[dyq])H*- 


(5.9) 


(5.10) 


(5.11) 


(5.12) 


As (5 —)• 0, A tends to Aq which is defined to be 


Ao 


Cc + 1 

2(ee-l)' 


The formula p5.12p shows that the resonance p5.11l) occurs only when Aq £ [—b, b] C (—1/2,1/2). 
This condition can be fulfilled only if Cc < 0. So, resonance occurs only when the dielectric 
constant of the inclusion is negative. A material with a negative dielectric constant is called a 
plasmonic material. 

Let 

Hq{t) := {d^q,£t[d,yq])n* (5.13) 

for ease of notation {fig is called the spectral measure). Since dfig is absolutely continuous, 
fi'q G L^{—b,b) (Lemma 14.3p . Observe that A = Aq + ic6 + 0{5‘^) for some positive constant c. 
So we assume from now on that 

A = Ao + i(5 (5.14) 

for simplicity and Aq G [—b,b]. Then (j5.12p can be expressed as 

fb u'(t) 
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In what follows we investigate how fast ||(/?(5||'K* tends to oo as <5 —)• 0 in a quantitatively precise 
manner. For that, we emphasize that 


1 (5 

vr (Ao - ty + (5^ 

is the Poisson kernel on the half space, and it holds that 

for all integrable function / if the limit exists. Here /(Aq—) indicates the limit from the left and 
/(Ao+) from the right. 


5.1 Generalized eigenfunctions 

For s > 0, let 5s be the Dirac mass at s and dehne 


dl : = 


V2pi{s) ^/"^Ss 
0 


and d::= 


y/2p2{s) ^^‘^{5s-5o) 


(5.17) 


Since K* is a multiplication operator given by (13.351) . d^ and are generalized eigenfunctions 
of K* corresponding to the continuous spectrum 77 ( 5 ) G cr(]K*), and and d?_^ are those 
corresponding to We emphasize that Pj{s)~^^‘^ appears in the definition of dj because of 

weights in the definition of W* and do is subtracted in the definition of d^ because /2(0) = 0 (if 
f G W* is continuous). Note that 


(f,d^)ty. = ^Pj{sy^^fj{s), j = 1,2. 


(5.18) 


Let 


^Pi:=U-yd^, j = l,2. (5.19) 

These are generalized eigenfunctions of fCg^. We emphasize that since d* ^ W*, ipi ^ T-L*. In 
fact, as we show below, \ijjl{x)\ tends to 00 and ipiix) oscillates fast as x approaches to one of 
the intersecting points. But, since 6s G 77“^ (M), we infer that ipi G Since Sq^i maps 

H~^{dQ,) into L^(dD), ((/j,= — ((/?,is well dehned for p G L^(dD), and we have 


for all ip G L'^idn). 

One can easily see that 

1 


(5.20) 


T-^A-yd^sm = 


20F 


Pi{s) 


So, we have from (13.161) 

ylix) = < 

Likewise, one can see that 
1.2 


2y vr 

ly/TT 


X G r+, 
X G r_. 


yl{x) = P 2 {s) ^/^/i(^'i(x),6'o)[e*^'^ii"^i - 1 ], xGF. 
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Since 'I'i(x) —)• oo as a: approaches to intersecting points of two circles, we can see from these 
formula that \'4>i{x)\ tends to oo and oscillates fast near the intersecting points. 

We can also compute the single layer potentials of these function. In view of ()3.24l) and 
(j3.25p we have 

Sanm = [ATMhCr^SUi^Pi] = -{ATMhC)-^[P{s)d^, 


and hence 


and 


Saniil^lKx) = < 


2y TT 
2^y TT 


X G r+, 
X G r_, 


SanlTp^Jix) = P 2 {s)^^^ _ i] ^ a; g r. 


Lemma 5.1 It holds that 


•Sani'^Plliz) 


1 


20F 

1 


P2{s)^^‘^ 


sinhs^2(z) 

sinh sOq ’ 

' C0shs^2(^) ^is^Uz) 
cosh sOq 



(5.21) 


(5.22) 


(5.23) 

(5.24) 


for z G \ n. 


Proof. Since ^I'l(.^) + i'i’ 2 (z) is analytic, the functions on the righthand sides are harmonic in 
M^\n. We see from (j5.21l) and (j5.22p that equalities in (|5.23p and p5.24p hold on F = dfl. Since 
'I'i(z) Pi^ 2 (z) ^ 0 as \z\ —^ oo, we see that functions on the righthand side of (j5.23l) and (|5.24p 
tend to 0 as \z\ —)• oo. On the other hand, since {tpi, 1 )-h* = Oj we infer that San['<pi]{z) —)• 0 as 
\z\ —)■ oo. So, (I5.23P and (15.241) for z G \ follow from the maximum principle. □ 

Let pq be the spectral measure defined by (j5.13p . Then we have from Lemma 14.31 that 


Pq{t) 


1 


|r/'(s)| 

1 

. h'(s)l 


\{duq,'ips)n*f + \{duq,ip-s)n*\‘^ , 


t G (0, b], 
t G [—6, 0), 


where s = rj ^(t) if t > 0 and s = rj ^{—t) if f < 0. We see from Lemma l5.ll that = 

Sani'f’s], and hence \{d,yq,'if^_g)-p* \ = \{d^q,ipi)n*\‘^ for j = 1)2. Thus it holds that 

25-1 (t) 

(5.25) 


Pqit) = 


W{rj-^it))y 

2ff2(-^) 


t G (0,6], 
t G [-6,0). 


Here we set for ease of notation 


9j{t) ■= 


, j = l,2. 


(5.26) 


Note that d^q belongs to is smooth except at intersecting points of two circles, real 

valued, and satisfies 

/ duqda = / Aq{x)dx = 0. 

Jon Jn 

So, gj is continuous on [—6,6] except at 0, and so is p'. 
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5.2 Plasmonic resonance 


In this subsection we show resonance occurs at the continuous spectrum in a quantitatively pre¬ 


cise way. Since fig is continuous in (— 6 ,0) U (0, b), the following proposition follows immediately 
from (15.161) . 


Proposition 5.2 Suppose Aq G {—b, 0) U (0, b). Let ips be the solution to (15.6p with A = Aq -|- i6. 
Then it holds that 


0—>-(J 


27r5ri(Ao) 

|r/'(^-nAo))| 

2vr52(Ao) 

|r/'(^-n-Ao))| 


if Xo G (0,6), 
if Xo G (-6,0). 


(5.27) 


Let us consider the case when Aq = 6 . We see from (15.151) . (j5.16p and ^5.251) that 


YimJWpsWn* 

()—>-(J 


TT 


hn{b-) = lim 




t^b- \v'{r]-^{t))\ 


Using P3.36I) we can see that 




(5.28) 


for some nonzero c. So we infer that 


lim 5\\ps\\n* = oo 
5^0 


as long as gi ( 6 ) ^ 0. We may rehne this result in the following way. 


Proposition 5.3 Suppose that Aq = —6 or 6. Let be the solution to (|5.6p with X 
Then it holds that 


Yimd^/^WpbWn. 


Cgi{b) ^fXo = b, 

Cg 2 {-b) */Ao = -6 


for some nonzero constant C. 


Aq + id- 
(5.29) 


Proof. Suppose that Aq = 6 . We see from (I5.15P that 






Since the integral over [—6,0] tend to 0 as 6 —)• 0, it suffices to consider the one over [0,6]. We 
have from (I5.25P that 


Jo (6-t)2-K62 Jo {b - t)'^ + \g'{g ^(t))] 

We show that there is a constant C 7 ^ 0 such that 

}^oJo (6-t)2 + 62 1r7'(r?-i(t))l^^“^- 
In fact, one can see from (|3.36p that 

1 _ 1 
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where e{t) is a continuous function on [0,6] and c is the number appearing in (|5.28l) . So, we 
have 

^o/o ib-tr + 5^wiv-Hm sTol (b-ty + 6^cVb^t 

By making change of variables, b — t = 5s, we have 


5^/^ tt , vr /■°° 1 

™o I (6 - t)2 + 62 “ c io (s 2 + l)Vi 


Using the constant C in (I5.30p . we write 
63/2 


I 


0 (6-t)2 + 62 


dt - Cg{b) 


I 


63/2 7r[gi(t)-ffi(6)] 

^0 (6-t)2 + 62 |7?'(r/-i(/))l 


dt + 5 ( 6 ) 


f 


63/2 


TT 


(6 - t )2 + 62 


dt-C7 


Using continuity of g and (I5.3U|) . it can be proved by standard estimates that terms on the 
righthand side above converge to 0 as 6 —)• 0. So (15.291) for Aq = 6 is proved. The case when 
Ao = — 6 can be handled similarly. □ 

Observe that the case when Aq = 0 is missing in Proposition’s of this subsection. We deal 
with this case when the source function is a polarized dipole in the next section. Results of this 
subsection show that if 

ffi(Ao)/0, (5.31) 

then resonance occurs at the order of 6 “^/^ or 6“3/'^^ namely Ilyp^H-^* —>■ 00 at the order of 6“^/2 
or 6“3/4 as 6 —>• 0. We will show in the next section that resonance occurs at Aq = 0 at the rate 
of 6 “" for any a < 1 where a varies depending on the location of the dipole source. We will 
also show that (I5.3ip is fulfilled with polarized dipole sources. 

It is worthwhile to mention that this order of resonance (at the continuous spectrum) is in 
contrast with that at the discrete spectrum where the order is 6 “^ (see a)- The following lemma 
shows that at the absolutely continuous spectrum (on any domain with Lipschitz boundary) 
resonance never occur at the order of 6 “^. 


Proposition 5.4 Let Q be a bounded domain with the Lipschitz boundary. Suppose that the 
NP operator on Ll has only absolutely continuous spectrum. Let (ps be the solution to (|5.6p with 
A = Ao + i5. Then, it holds that 

Iim6||<y95||^* = 0. (5.32) 


Proof. Let g be the spectral measure (as defined in (|5.13l) i for the NP operator on 11. Since 
only absolutely continuous spectrum is nonempty, g' € L^(M^), and it holds that 



h'jt) 

(Ao - t)2 + 62 


dt. 


(5.33) 


For simplicity, we assume Ao = 0. We decompose 62 ||(/? 5 |||,* as 


f°° 62 r 

' —CXD ^ “ 1 “ ^ 


00 « 


(5<|t|<2''5 


-^0 + Ik- 


k=l 


Since g' is integrable, we have Jo 0 as 6 —)• 0. For k > 0, we have 


h = 


f 


^g (t) dt < 


1 


-I5<|t|<2'=(5 1 + 22(*' ^W2''-i(5<|t|<2''5 


g'{t) dt. 
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Observe that 


Ij!{ t) dt ^ 0 

/2'=-i(5<|t|<2*:5 

as (5 —)■ 0. So by Lebesgue dominated convergence theorem, we have 


/ 

J 2 '= 


E 

k=l 


1 


1 22 (*^- 1 ) 


/ 

J 2 ''-l 


5<\t\<2'^S 


(it —)• 0 


as (5 —)• 0. This completes the proof. □ 


6 Resonance by dipole sources 


Suppose that the source / is given by 

f{x) = a-V5^{x) (6.1) 

where a is a constant vector and z is a point outside 11. In this case, we have 

q{x) = a ■ V 2 ;r(x — z) = —a ■ V 2 r(x — z). (6-2) 


So we have 

{duq,^i)n* = [ -^T{z - y)Sdn[ipi]{y)da{y). (6.3) 

Jdn 

Note that 

/ ^r (2 - y)San[i^l]{y) da{y) = VQnS^n[^|Jl]{z) 

Jan (J^y 

Recall that the double layer potential T)qq, satisfies 

^9oMI+ = (“ 2 ^ + ^an)W\ 

for all (f € L‘^{dVL). So, we have from (j5.20l) 

{‘f,TJ'dnSan[i>i]\+) = A ~ 2^ + l^dn)Sdn['^Pi]) 

= (( “ 2 ^ + ^dn)['f]^'Sdn['4’i]) 

= - >^dn)[‘P],i^i)'n* 

= ( - ^ + (-l)^'^S(s))(i?,>5ao[V'i]). 

So we have 

'UdnSani'ipi] = ( “ ^ + i-'^y^^v{s))San[ipi] on dn. 

Since 'DQQSQQ['ipi]{z) and Sanliplyz) are harmonic in \ 11 and decay to 0 as \z\ 00 , we infer 

that 

[ -^'^{z-y)San['4’i]{y)da{y) = (- ^ + (-l)^+^T/(s))5ao[V’^](z). 

Jan (^’^y ^ 
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We then see from (| 6 . 3 p that 


9j{vis)) = = 

Choose z € \ 0 so that 




|a-VSan[V>i](2)r 


^ 2 ( 2 :) / 0 , 

which is equivalent to that z is not on the real axis. We see from Lemma l 5 .II that 

(a • V'I'2(2))^ cosh^ 'S'I'2(^) + (a • V'I'i(z))^ sinh^ s'I'2(2) 


1 sV('S) 


47r sinh^ sOq - 
Since 'I'lPz) + i'^2iz) is conformal, we have 

(a-V 4 'i(z ))2 + (a-V^2(2))V0 


( 6 . 4 ) 

( 6 . 5 ) 


( 6 . 6 ) 


for any a. We choose a so that 

a-V'Li(2;)/0 and a • V'L2(2:) / 0 . ( 6 . 7 ) 

Then, |a • V 5 an[V’s](- 2 )l / 0 foi' s, and hence 5 i(t/(s)) / 0 , or (I 5 . 3 ip is satisfied. Similarly, we 
obtain 


|a-V5aQ[V’"](2)|" 
1 S^P2is) 


dvr cosh^ sOq 


(a • V'I>2{z)f sinm s'I'2(2;) + (o • VTi(2;))^ cosm s'I'2(^) 


and hence g2{'i]{s)) / 0. 

Let us now consider the resonance at Aq = 0 . To do that we need to look into the behavior 
of gj{t)/\r]'{r]~^{t))\ as t ^ 0 +. We see from (j 6 . 4 p and ( 16 . 6 p that 

(5(5)) = (^^se2.(-0o+|f 201) [1 + Q(g-2|vi/20h)j as s ^ 00 
for some constant Ci 7^ 0 . We also see from the dehnition (| 3 . 31 l) of g that 


g{s) = C2e-2^°"[l + 0(e-2(^-2'^o)")] 


and 


r/'(s) = -200^26-2®°® [1 + O(e- 2 (^- 20 o)")] 
as s —>■ 00 for some constant C2 7^ 0 . Therefore we have 

for some constant C3 7^ 0 and e > 0, which implies that 

__= C'il + e^ (t) 

as t —)• 0+, where the error term ei(t) satisfies 

|ei(t)| < C|logt|t-l^2OI/0o+., 


( 6 . 8 ) 


( 6 . 9 ) 
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for some constant C. Likewise we obtain 


__= Cd\ loetlt“l'^2L)l/0o _|_ Q (f) 

as t ^ 0 + for some / 0 , where 62 (t) satisfies 

|e 2 (t)l < C\ (6.10) 

Theorem 6.1 Suppose that the source f is given by (EH) with a and z chosen to satisfy (EH 
and ()6.7|1 . Let (ps be the solution to ()5.6p with \ = id. Then there is a constant C > 0 such that 

I log(i|-i(ii+l^2(.)|/eo||<^^||^. ^ C (6.11) 


as 5 —)• 0 


Proof. We have from (15.121) 


/. 

rb 1 

Jo + 


W l_^ ^2 + 




/o t'^ + 6‘^2\r]'{r] i(t))| 

=:h{6) + l2{6). 


dt “h 


T^92it) 


lo t"^+ 6^2\r]'{p i(t))| 


dt 


We show below that 

I log5ri(5^+l'^2(z)|/eoj^.(^) ^Cj^O (6.12) 

as (i —)• 0 for j = 1,2. Then (16.111) follows immediately. 

Let us prove ()6.12l) for j = 1. The case for j = 2 can be handled with in the same way. Let 
Cl be the function appeared in ()6.8p . Then one can see from ()6.9I) . by the scaling t = 5s, that 


lim |log5r^5^+l'^2(^)|/eo f 

^^0 Jq 


ei{t) 


t"^ + (5^ 

On the other hand we see by the same scaling that 

Jo t + 5 Jq 


dt = 0. 


(6.13) 


s^ + 1 
<5 rb/S 


log s 


+ 1 


log 5 

rd 

/ +/ =-.Ih{d) + Il2{5). 

Jo Js 
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If s > 5, then log s/ log d —)■ 0 as 5 —)• 0 for all s. So by Lebesgue dominated convergence theorem 
we have 

poo 1 

1/2(5) ^ - g-\^ 2 {z)\/eo ^ Q ( 6 . 14 ) 

Jo + 1 

as 5 —)• 0. If 0 < s < 5, then | log s/ log 5 + I| < 2| log s\ if 5 is sufficiently small, and hence 

Ih[5) <2 [ logs|s-l'^2(^)|/0o 

Jo ■s^ + l 

Since |'I' 2 ( 2 )| < ^Oj the integrand on the righthand side is integrable, and it follows that 

//i(5)^0 (6.15) 
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as 0. By combining (|6.14p and (IG.lSp we infer that 


Jo 
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+ 6 “^ 


logt|fl'^2(2)|/0o^^ ^ 


Jo + 1 


as 0, which together with (I6.13P implies (I6.12p for j = 1- This completes the proof. □ 
Estimate ph.lll) shows that resonance occurs at Aq = 0 at the rate of | log ( 5 |^/ 2 ( 5 -i/ 2 -|'i' 2 (^)|/ 20 o_ 
So, the stronger resonance occurs if \^2{z)\/9o is closer to 1, in other words, if z is closer to 
517. But since |T 2 ( 2 :)| < 9q, the resonance at the rate of is never achieved as we showed in 
Proposition 15.41 
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